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A multi-mode model of a non-classical atom laser produced by outcoupling from a 

Bose-Einstein condensate with squeezed light 

S.A. Haine and J.J. Hope 

Australian Centre for Quantum-Atom Optics, Department of Physics, 
Australian National University, Canberra ACT 0200, Australia 

We examine the properties of an atom laser produced by outcoupling from a Bose-Einstein con- 
densate with squeezed light. We introduce a method which allows us to model the full multimode 
dynamics of the squeezed optical field and the outcoupled atoms. We show that for experimentally 
reasonable parameters that the quantum statistics of the optical field are almost completely trans- 
ferred to the outcoupled atoms, and investigate the robustness to the coupling strength and the 
two-photon detuning. 

PACS numbers: 03.75.Pp, 03.70.+k, 42.50.-p 



I. INTRODUCTION 

Certain precision measurements are improved by using 
slow- moving massive particles 0. In a Sagnac interfer- 
ometer, for example, the inherent sensitivity of a matter- 
wave gyroscope exceeds that of a photon gyoscope with 
the same particle flux and area by 11 orders of mag- 
nitude 0. The high spectral flux and associated first 
order coherence of atom lasers make them an obvious 
choice for the input of such devices. Although current 
atom laser experiments usually operate in a regime lim- 
ited by technical noise, the fundamental limit of these 
measurements will be caused by the shot noise of the 
atomic field, which will be intrinsic to all interferome- 
ters without a non-classical atomic source. Sensitivity 
is increased in optical interferometry by 'squeezing' the 
quantum state of the optical field, where the quantum 
fluctuations in one quadrature are reduced compared to 
a coherent state, while the fluctuations in the conjugate 
quadrature are increased. In the context of atom optics, 
it is therefore interesting to ask whether highly squeezed 
atom optical sources can be produced. 

Generation of squeezed atomic beams has been pro- 
posed by either utilising the nonlinear atomic interactions 
to create correlated pairs of atoms via either molecular 
down conversions or spin exchange coUisions Q, Q, 0, 
or by transferring the quantum state of a squeezed optical 
field to the atomic beam , , |1| , . 

The generation of nonclassical light is well established 
experimentally This suggests that a nonclassical 

atom laser output could be generated by transferring the 
quantum state of an optical mode to an atomic beam. 
Moore et al. showed that a quantized probe field could be 
partially transferred to the momentum 'side modes' of a 
condensate consisting of three-level atoms in the presence 
of a strong pump field, via a Raman transition [g. Jing 
et al. performed a single mode analysis of the atom laser 
outcoupling process for a two-level atom interacting with 
a quantized light field, and showed that the squeezing 
would oscillate between the light field and the atomic field 
at the Rabi frequency 0. Fleischhauer et al. IH showed 



that Raman adiabatic transfer can be used to transfer 
the quantum statistics of a propagating light field to a 
continuously propagating beam of atoms by creating a 
polariton with a spatially dependent mixing angle, such 
that the output contained the state of the probe beam. 

It has been demonstrated that the complicated 
multimode dynamics in atom laser outcoupling can cause 
effects such as back coupling and a 'bound state' signif- 
icantly limiting the flux. Although it is well established 
that a system of three level atoms interacting with a 
quantized probe field via a Raman transition can exhibit 
some degree of quantum state transfer, what remains to 
be demonstrated is whether the multimode effects inher- 
ent in atom laser outcoupling will inhibit the clean trans- 
fer of the quantum state to the atom laser beam. In a 
previous paper we showed that even when we assume 
that both the probe field and condensate mode are single 
mode, the outcoupled atoms can still exhibit complicated 
multimode dynamics. For complete quantum state trans- 
fer, the the atoms will have to undergo a quarter Rabi 
oscillation in the time taken to leave the coupling region, 
and the light from the probe field will have to be com- 
pletely absorbed. In our previous paper we assumed a 
single mode probe field, which is a valid approximation 
if the light makes many passes through the condensate 
before it is completely absorbed, as in a high finesse cav- 
ity. For a more realistic single pass experiment, a multi- 
mode model of the light field is required. In this paper 
we develop a technique which allows us to model the full 
multi mode dynamics of the probe beam and the outcou- 
pled atoms, and allows us to investigate how multimode 
dynamics affect the quantum state transfer. 



II. MODEL 

We model an atom laser in one dimension as a BEC 
of three-level atoms coupled to free space via a Raman 
transition, as shown in Fig. ^ The optical field affect- 
ing the 1 2) 1 3) transition (pump field) is assumed to 
be strong and is well approximated by a monochromatic 
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FIG. 1: Internal energy levels of a three level atom. A con- 
densate of state |1) atoms confined in a trapping potential is 
coupled to free space via a Raman transition. The two fields 
of the Raman transition are a probe beam (annihilation oper- 
ator E{x,t)) and a semiclassical control field {Q,2s,{x,t)) that 
is detuned from the excited state by an amount A. 



classical field ^23 (a;, t) = f]23e'^'='"'"^"°"^^*^ where rJss is 
the single photon Rabi frequency and A is the detuning 
from the excited state. 

The Hamiltonian for the system is 



'H — Ti-atom + Ti-int + "Hlight 



(1) 



r fc2 
J 2m 

J 2m 

+ fi- j {'^2{x)ipl{x)Vl{x,t) +h.c.)dx 

+ %13 / {E{x)'ii}i{x)'4>l + h.c.)dx + Uught 



where Vl2^{x,t) = ^^^(,^{ko^-{'^o-A)t) ^^^^^ 
rabi frequency for the |2) — > |3) transition, Hq is the sin- 
gle particle Hamiltonian for the trapped atoms, and m is 
the mass of the atoms, i/^iix) is the annihilation operator 
for the condensate mode (internal state |1)), il^six) is the 
annihilation operator for the atoms in the excited atomic 
state (1 3)), and ip2{x) is the annihilation operator for the 
untrapped, freely propagating atoms (|2)). The atomic 
field operators satisfy the usual bosonic commutation re- 
lations: 

IMx), M^')] = [4>Ux), ViI(x')] = o, 



£^(0;) is the spatially dependent annihilation operator for 
the probe field, satisfying 



[Eix), E^x')] ^ d{x - x') 



(3) 



The coupling coefficient, 513 = ^y^^, where ^13 is 

the dipole moment for the |1) — > |3) transition. We have 
assumed that 313 (wfc) is approximately flat in the range 
of interest of our system. 

The equations of motion for the Heisenberg operators 
are: 



iipiix) 
#3(2;) 



-^■ipi{x) + gi3ijj3{x)E^x) 



2m 



2m 

gi3E{x)ipi{x) 
d 



6ijS{x ■ 



(2) 



iE{x) = - {ujQ - A))E{x) + gi3ipi(x)ipl{x) 

(4) 

where i^six) = i^ae^^"""^)* and E{x) = ES'^o-^^)t^ ^he 
population of the excited state |3) will be much less 
than the other levels when the detunings of the opti- 
cal fields are larger than the other terms in the sys- 
tem (including the kinetic energy of the excited state). 
Furthermore, most of the dynamics will occur on time- 
scales less than -g-, so in this regime we can set ^^3(2;) = 
_I^gifcoa;^2(a:) - ^E{x)'ipi- If the condensate is ap- 
proximated by a coherent state containing a large number 
of atoms and the outcoupling is weak (ie. the condensate 
wave function is not strongly perturbed by the outcou- 
pling) we can write ?Ai(a;,t) = •\/iV0o(2^)e~*"'*, where 
4>q{x) is the ground state wave function of a harmonic 
trap with trapping frequency Ut , and N is the mean num- 
ber of atoms in the condensate. The assumption of ignor- 
ing the backaction of the outcoupling on the condensate 
wavefunction is only valid if the number of photons in the 
probe beam is much less than the number of atoms in the 
condensate, as it will need to be in a successful quantum 
state transfer experiment. We have ignored the atom- 
atom interactions in our model, which will only be valid 
if the condensate is dilute. Strong atom-atom interaction 
will cause complicated dynamics of the quantum state of 
the condensate mode, which may not necessarily inhibit 
the efficient generation of a squeezed atom laser, but we 
can not model this evolution with the techniques used in 
this paper. Inclusion of such effects would require a more 
complicated technique, such as a phase space method 
|l2j |. With these approximations, our equations of mo- 
tion for the outcoupled atomic field and the probe field 
become 

ii:{x) = Hai}{x)^ilc{x)e~'^°''E{x) (5) 
iE{x) = HbE{x)-Q.*c{x)e''">''i>{x) (6) 
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with H, 



2m dx'^ 



" dx 



^100 Wl^ - - A)), nc{x) = S^ViV0o(x), and 

l/, = '026*'^**. 

Equations of the form of (jSJ and are common in 
quantum atom optics problems. In the next section we 
introduce our technique for solving these equations. 



III. SOLUTION METHOD 

At t = the field operators can be expanded as 

i 

E{x,t = Q) = "^p,{x)bi (8) 

i 

where fi{x) and Pi{x) represent an expansion in any 
orthonormal basis, the operators and bi represent 
Schrodinger picture annihilation operators for the ith 
mode of the atomic and optical fields respectively, and 
a|ai (blbi) represent the number of atoms (photons) in 
mode i. From this we can postulate that the solution to 
Eq. © and Eq. © is 



'4'ix,t) = ^fi{x,t)ai+^gi{x,t)bi 



(9) 



E{x,t) = '^pi{x,t)bi +'^qi{x,t)ai (10) 

i i 

By substituting this ansatz into Eq. jSJ) and Eq. ^ wc 
obtain equations of motion for the mode functions fi{x), 
gi{x), pi{x) and qi{x) is given by 

iUx) = Hah{x)-n,e-'^'>-q,{x) (11) 
ig,{x) = Hagr{x)~^,e~'''°^p,{x) (12) 

iMx) - Hw^{x)~Kf^'''"''g^{x) (13) 

lUx) = Hbq,{x)-n:e'''^"'f,{x) (14) 

with gi{x,t = 0) = qi{x,t = 0) = 0. In practice we 
are free to choose any initial conditions we like for the 
fi{x)s and pi{x)s, as long as they form an orthonormal 
basis. From the solutions to these equations we can ob- 
tain the dynamics of any observable of the system. How- 
ever, keeping track of all the mode functions (particularly 
when a numerical solution is required) can be daunting. 
In our system we initially have no outcoupled atoms, 
and if we choose our basis carefully, we can choose it 
so that the photons only initially occupy one mode. i.e. 
|*(t = 0)) = \0,0,..0) atoms® \l,0,...,0)Hght whcrc I7) 
represents an arbitrary state of a single optical mode. In 
our subsequent calculations we choose the mode basis to 
be plane waves, with bo representing the annihilation op- 
erator for a plane wave with momentum kp. This means 
that the initial condition of ^0(2;) is a plane wave with 
momentum k„. 



By noticing that a^, for all i and bi, i ^0 acting on 
our state return zero for all time, we can write Eq. 
and Eq. H1U|I in a more illuminating way. 

-0(a;,i) go{x,t)bQ + V^{x,t) (15) 
Eix,t) = poix,t)bo + VEix,t) (16) 



with 



V^ix, t) 



'fi{x,t)ai+^gi{x,t)bi (17) 



VEix,t) = ^qi{x,t)ai + ^pi{x,t)bi (18) 



Instead of solving for the coefficients of the annihila- 
tion operators in and Ve individually, we note that 
since we know how they operate on our state for all time, 
the only extra information we need are their commuta- 
tion relations. By noticing that [ilj{x,t), ip''{x',t)] = 
[E{x,t), E'^{x',t)] = 5{x - x') and [a„ a]] ^ 

[bi, foj] — 5ij with all other commutators zero, it is easy 
to show that 

[V^{x,t),vl{x',t)]=5{x-x')-gl{x',t)go{x,t) (19) 



[VE{x,t),vl{x',t)]=5{x-x')~pl{x',t)po{x,t) (20) 

From the solutions for go{x, t), po{x, t) and these com- 
mutation relations, we can calculate any observable of 
the system. This method converts the full multi-mode 
quantum field system into two coupled classical modes 
that change dynamically with time. 

To reduce computational burden when solving for 
go{x,t) and po{x,t) numerically, we have transformed to 
the momentum shifted variables po (2;) = po{x)e~'^^^''^ and 
gQ{x) = gQ{x)e^''^°~^'p'^ giving us 

= (-:^V2- A|ko-kp|V 
Zra m 



2m 



|ko 



'231 



A 



(21) 

ujt)go{x) - ^cPo{x) 



ipo{x) — (— icV + 5 — 



I5l3| 



■N\Mx)\^)Po{x) 



ni{x)go{x) 



with S 



c|kp| - c|ko| (since c|ko| 



LOq 



(22) 
A) be- 



ing the two photon detuning. Physically, kp repre- 
sents the momentum of the initial mode of the opti- 
cal state with non-zero occupation, as the initial con- 
dition on po{x) defines the mode that is occupided. Op- 
timum coupling between the atomic and optical fields 
will occur when they are resonant. This will occur when 

S^^\ko- kpp + ^|0o(O)P -\^-u;,^6o. 

By noticing that the evolution of the optical field is 
trivial at high frequencies (of order c|kp|), and the evolu- 
tion we are interested in (the interaction with the atomic 
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field) will occur on frequency scales of less than ftc, we 
can make the slow envelope approximation and ap- 
proximate the dynamics of po{x,t) as 



= {6~^-^\M^r)Po{x) ~ nigoix) (23) 

In the following section we will look at the solutions to 
these equations and use them to look at some properties 
of the outcoupled atoms. 



IV. PROPERTIES OF THE OUTCOUPLED 
ATOMS 

The density of outcoupled atoms and the optical den- 
sity (mean number of photons per unit length) are given 

by 

{EHx)E{x)) = |p(x)p(6t6„) 

The observables {ij:'^ {x)'il}{x)) and {E'^ {x)E{x)) are only 
related to our solutions for g^ix, t) and Pq{x, t) by such a 
simple expression because of the specific choice of the 
basis for the optical field. We solved Eq. (|21|) and 
Eq. H23() numerically for goix^t) and po{x,t) using a 
4th order Runge Kutta algorithm with a cross propa- 
gation step using the numerical package XMDS jlJI . We 
chose parameters realistic to atom optics experiments 
with Rb*^ atoms. Unless stated otherwise, we have set 
TO = 1.4 X 10-25 kg, LJt = 20.0 rad s-\ gig = 28.9 rad 
s-^, N = 10^ and A = 10" rad s"^ Fig. Elshows typical 
results for the densities |go(a^)P and |po(2;)P after t — 7.2 
ms of evolution. 

For optimum transfer of squeezing from the optical 
field to the atomic field, the quantum efficiency of the 
outcoupling process (i.e. number of atoms per photon in 
the probe beam) will have to approach one. When this 
is the case, the probe field will be completely absorbed. 
Fig. 121 shows that overcoupling can significantly reduce 
this efficiency. 

Reduction in the variance of the flux of an atom laser 
beam would be a measure of how much squeezing is trans- 
ferred to the atom laser from the optical field. However, 
in a fully multimode model, the variance of the flux is 
infinite. Instead, we use the density integrated over a 
small region as a measure of how 'quiet' the atom laser 
beam is. We define the operator 



N : 



ijj^ {x)iij{x)dx^ 



(24) 



which represents the number of atoms in a region of space 
between xi and X2 in the path of the atom laser beam. 
Physically this operator would be a measure of the 'in- 
tensity' of the atom laser beam. The variance of N is 
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FIG. 2: \go{x)\'^ (dashed line) and |po(a;)P (solid line) at 
t = 7.2 ms as found numerically for (a) S — So and 5123 = 
2.1 X 10^^ rad s-\ (b) 5 = So and ^23 = 3.2 x 10^^ rad 
and (c), S = So + 4 x rad s"^ and ^23 = 2.0 x 10^^ 
rad with all other parameters given in the text. The 
densities are related to the functions \go{x,t)\'^ and |po(a;,t)P 
by a factor of (bgho), ie. the initial optical density of the 
field. The dotted line (|(^o(2;)|^) represents the density profile 
of the condensate (density in arbitrary units). In (a), the 
probe beam is attenuated as it passes through the condensate, 
producing a steady flux of atoms. In (b) the coupling strength 
is too strong for ideal quantum transfer, causing significant 
back coupling of atoms into the condensate, which causes a 
bound state of atoms and reduces the flux. In (c), the two 
photon detuning is such that the light in the probe beam is off 
resonant with the Raman transition. An initial pulse of atoms 
is ejected from the condensate, then the flux reduces as atoms 
remain bound to the condensate, poix) was normalized at 

0.02 



t — such that 



-kol 



J^i^ \p{x)\ dx — 1 with 



V{N) = (N^) - (iV)2 

{-ip^ {x')il;{x')'ip'^ {x)il>{x))dxdx' 

{'4jHx)i;{x))dxY 

{'ip^x')^'^ {x)ij{x')i;{x))dxdx' 

+ / {ii;\x)^j{x))dx -( I {ip\x)ip{x))dx^ ^ 



Using our solution for ^{x,t) this becomes 



V{N) = Nl(^{hlhlbok) ~ {blho?) + Ngiblbo) (25) 

where Ng = JJ' \g{x)\'^dx. To compare this to the ini- 
tial amount of squeezing in the optical field, we define 
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the operator Nq = J^,^ E\x)E(x)dx, where x'2 — x'^ — 
ri|k™-ko| ~ ^1)1 is. the range of integration for the 
optical fields is larger than for the atomic field by a fac- 
tor of the ratio of the light speed to the mean atomic 
speed. This is because if we have a well defined number 
of photons in a region of length L initially, if there is 
ideal quantum transfer, the number of particles will now 
be confined to a region of length ^\^p~^o\ T\ie variance 
at t = is then 

V{Na) = N^(^{blblbobo) - {blk)^) + Npiblk) (26) 

where Np = \p{x,t = 0)\^dx = 1, since we have nor- 
malised Pq{x, t — 0) such that represents the num- 
ber of particles per length X2—X1 . Since we have assumed 
that the initial state of the occupied mode for the optical 
field is a plane wave, the variance is homogenous along 
the length of the beam at t = 0. The normalized variance 

viNo) = ^ is then 



viNo) 



{{blblbobo) - (blbo) 



{blbo) 



(27) 
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FIG. 3: vpock versus time for (a) S — So and fl23 = 2.1 x 10^^ 
rad s~\ (b) S = So and ^23 = 3.2 x 10^^ rad s'\ and (c) 
5 = (5o + 4 X 10^ rad s'^ and ^23 = 2 x 10^^ rad s^^ As 
the flux of the atom laser beam becomes steady, the vpock 
reaches it's minimum value, which is limited by the quantum 
efficiency of the outcoupling process. The detector region was 
chosen such that xi — 0.04 mm, X2 = 0.06 mm. 



In terms of v{No), v{N) 



IS 



viN) = NgviNo) + (1 - Ng) 



(28) 



As Ng 1, v{N) v{No). If the initial state of the 
optical field had perfectly well defined number (i.e. a 
Fock state), when?;(7Vo) — 0, andw(7V) = l—Ng = vpock- 
If a coherent (or vacuum) state is used to outcouple then 
v{No) = 1, regardless of the efficiency of the outcoupfing 
process. Fig. j3] shows vpock versus time for the cases 
shown in Fig. 

To get complete quantum state transfer (and hence the 
minimum possible variance in the flux) the coupling be- 
tween the optical and atomic field will have to be strong 
enough such that the optical field is completely absorbed. 
However, if the coupling is two strong, there will be sig- 
nificant back coupling of the atoms into the condensate 
as was seem in reference [IH, We can estimate the 
optimum coupling desired by equating the average time 
taken for an atom to leave the coupling region (conden- 

-m/{h(kp 



sate) Ti 



ko)) where 



- , , v.ii^.±v. ^ , — ^ is the 
spatial width of the condensate, with the quarter period 
Rabi oscillation Sapii = „ ^ r,^/ -.^ ■ This will occur when 
^23 ~ 2.3 x 10^2 j.a,d s-i. 

Obviously, the amount of squeezing that can be trans- 
ferred to the atom laser beam is strongly dependent on 
the strength of the coupling and the two-photon de- 
tuning. These factors introduce complicated multi-mode 
behaviour into the atom laser beam which reduces the 
fiux and also will broaden the line width, as the atomic 
mode becomes less monochromatic. Fig. 0] shows the 
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FIG. 4: (a) The minimum value of vpock as a function of 5— So 
for $123 ~ 2.15 X 10^^ rad with all other parameters given 
in the text. Optimum squeezing occurs at around S — So = 800 
rad s~^. (b) Minimum value of vpock as a function of 0,23 for 
S — So — 800 rad s^^. Optimum squeezing occurs at around 
0.23 = 2.2 rad s"\ 



minimum value of vpock obtained for different values of 
S and n23- 

Our numerical simulations show that at the appropri- 
ate coupling strength and two photon de-tuning, vpock 
tends to a very small number < 0.01, indicating that the 
squeezing from optical field is completely transferred to 
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the atom laser beam. However, we cannot give a quan- 
titative limit, as we are limited by round off error in our 
calculation of Ng. Obviously for complete transfer of the 
squeezing to occur, the light would have to be completely 
attenuated as it passes through the outcoupling region. 
Our simulations show that when = 2.2 x 10^^ rad 
s~^, then the squeezing is transferred with an efficiency 
greater than 0.99 while the optical field is attenuated by 
a factor of 10**. In any realistic experiment the squeez- 
ing measured in the atom laser beam would be limited 
by other factors such as limited squeezing in the optical 
beam and detection efficiencies, rather than by this level 
of quantum state transfer. 



V. CONCLUSION 

We have shown that a multimode model of atom laser 
out-coupling can display a very high level of quantum 
state transfer when outcoupled under appropriate con- 
ditions. This transfer is degraded when the system is 
overcoupled or off- resonant. The conditions necessary for 
optimum quantum state transfer are also the once which 
give the maximum clean flux for the atom laser beam. 

This research was supported by the Australian Re- 
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